Hikami boxes and the Sinai billiard. 
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Diagram, known in theory of the Anderson localization as the Hikami box, is computed for the 
Sinai billiard. This interference effect is mostly important for trajectories tangent to the opening 
of the billiard. This diagram is universal at low frequencies, because of the particle number conser- 
vation law. An independent parameter, which we call phase volume of diffraction, determines the 
corresponding frequency range. This result suggests that level statistics of a generic chaotic system 
is not universal. 



Any function, which describes the physical properties 
of a single particle bounded system, is called universal in 
this work, if it is a function of either frequency or time 
and has the only one parameter: the mean level spacing 
A divided by H. 



I. FORMAL DERIVATION. 

Following Ref. |l| the quantum evolution of the parti- 
cle inside the two-dimensional billiard is convenient to 
describe by the scattering matrix S(9,9'). For given en- 
ergy E = h 2 k 2 /(2m) the wave numbers of the incident 
and scattered wave functions are fcj = (k cos(0), k sin(0)) 
and kf = (kcos(0'),ksm(0')) correspondingly. The wave 
reflected n times from the billiard walls is described by 
the n-th power of the matrix S. 

In the particular case of the Sinai billiard one needs 
the product of two matrices S = S^So j where is the 
scattering matrix of the square, and is the scattering 
matrix of the circle. A similar decomposition was used 
for quantization of the Sinai billiard.H 

The direct product of two scattering matrices is a den- 
sity evolution operator. The two ways to pair arguments 
of two matrices give diffusion and cooperon-like opera- 
tors: 



= £< 
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where / = [r x k] z is the angular momentum measured in 
units of h, the z-axis points out of the plane. 

Operator D is a one step classical evolution operator. 
It is well defined if position of a classical particle on the 
energy shell of the phase space f, k/k is described by I, 0. 
This is the case for the Sinai billiard, see Fig. [I], if we 
agree to compute / and 6 for a particle, just before it hits 
the square. Thus each point in the phase space corre- 
sponds to either straight piece of trajectory, for example 



the point A, or trajectory with one reflection from the 
circle, for example the point B in Fig. [l| 

Each trajectory in the configuration space is a sequence 
of straight lines (segments), it becomes a sequence of 
points in the phase space; the reflection law generates 
a map in the phase space 



2 = M(1), l = {6 u h), 2 = {0 2 ,l 2 ) 



(2a) 



The density evolution operators Eq. ([!]) computed within 
semiclassicalH and diagonalB approximations are C ~ C, 
and T> w _D, where 

C£(l; 2) = D£(l; 2) = 27re Mtl ' 2 (5(2 - M n {l)) , (2b) 

A4 n means n iterations of the map M, and t\ 2 is the 
time of flight along the trajectory. 

The diagonal approximation is justified when classical 
actions of different trajectories are not correlated. The 
action of the trajectory is proportional to its length. The 
only scale oi the length-length correlation function is the 
system sizeQ, and the scale of the action correlations is, 
therefore, HNh, where Nh is the number of open chan- 
nels in the system, it is the effective dimensionality of the 
scattering matrixtl, and it is the analog of the Heisenberg 
time for maps. This argument provides us the condition 



n < N H or u > A/h , 



(3) 



where A is the mean level spacing. 

Under the condition Eq. (|^) the density-evolution oper- 
ator must preserve the invariant measured, and therefore 
J2 n tr^S must have single pole at uj = or alternatively 



1 



(4) 



This sum rule holds for the Frobenius - Perron operator! 
Eq. J2b|). In this case it isJsnown also as the Hannay - 
Ozorio de Almeida sum ruleo, and characterizes decay 
of Frobenius - Perron modes. 

Dashed lines in Fig. [l] mark the parts of the phase 
space where the semiclassical approximation fails because 
of diffraction. The phase volume of diffraction is much 
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smaller than the phase volume of the system. There- 
fore V w D + SD, and the first order correction to the 
evolution operator is 



« n— 1 
*D2(1;2)= / d3d4 53^'(l;3) 

n' — 1 



F D (3;4)-<5(3;4) 



(3; 2) 



F D (1;2) = 5 u J D (h 



fo(l,0)= / dl'S l+l , /2 S[_ ll/2 e 



Wo 



(5) 
(6) 
(7) 



where 5(1; 2) = 2 7 rS h iJ(9 1 -9 2 ), 3 = (0 3 , Z 3 ), 4 = (0 4) Z 4 ), 
Jd3 = Y. h I% S, = -H < t 1 \kR)/H < t 2 \kR) is the 
scattering matrix of a circle. The diffraction coefficient 
foil, 9) for kR-a' < I < kR + a", where a' ~ «" ~ 
(fci?) 1 / 3 , may be approximated in a number of wayso, see 
Fig. [|. Outside this interval, for Z < kR — a' , it must 
reproduce the map fo{l,9) — 2tt5(9 — arccos(^)), and 
for l> kR + a" one has f D (9) = 2ttS{9). 




FIG. 1. Configuration and phase space of the Sinai bil- 
liard. Four areas in the phase space correspond to the four 
sets of trajectories hitting four walls of the square. Reflec- 
tions from walls generate a map. It takes the particle from 
one area of the phase space and put it to other area. 

In order to compute the so-called interference correc- 
tion one should examine all possible pairing of arguments 
in the product of S matrices. The further computation 
is eventually the same as in Ref. || and the correction 
mixing two diffusion and one cooperon like operators is 



8D n (l,2) = j d3d4 (8) 
x |d(1; 3)£>(4; 2) [Fj(3; 4)C(4; 3) + F/(3; 4)C(3; 4) 
+ D(l; 3)£>(3; 2) [F D {3; 4)C(4; 3) + F D (4; 3)C(3; 4) 
+ D(l; 3)£>(1; 2) \f d (3; 4)C(4; 4) + F D {3; 4)C(3; 3) 



where , 3 = (9 3 + tt, -Z 3 ), 4 = (0 4 + n, -Z 4 ). The right 
hand side of Eq. (||) should be understood as a sum over 
all possible powers n\ + n 2 + n% = n of the density evo- 
lution operators, i.e. 5D n cx D ni C ,l2 D n \ The diffrac- 
tion kernels are defined in the vicinity of tangency to the 
cylindrical mirror: 



Fi(h,9i;l2,92) = 5iii 3 /i(ii, 9\ — 9 2 ) 

fii l ,8) = J dl'Si +v / 2 S-i + v/2e 



n'o 



(9) 
(10) 



where Si is the same as in Eq. (0), and 2Re/j for kR = 50 
is shown in Fig. || In order to avoid difficulties with self- 
tracing trajectories we put F/(l;2) = Fn(l;2) = for 
h <kR — a' and h > kR + a". 

The choice of the constants a', a" is restricted by the 
particle number conservation law 



d3dA F J (3;4) + F;(3;4) + F D (3;4) 







(11) 



Indeed, the interference correction Eq. (g) has a form of 
the interaction of the Frobenius - Perron modes. There- 
fore, the sum over all diffraction and interference correc- 
tions shifts the u> = pole of the evolution operator. In 
order to avoid the contradiction with the particle number 
conservation law in the form of Eq. (|^) the shift must be 
less or of the order of the mean level spacing A. (In 
this case the shift should be neglected because of the 
condition Eq . (pj ).) The shift of the pole depends on the 
integral Eq. (jllj). After all this integral must be so small, 
that we put it equal to zero. 



Id 



19 19 
FIG. 2. The numerical calculation of the diffraction and 
interference kernels. 

The interference correction Eq. (||) contains a small 
parameter: the probability for a given trajectory to visit 
the diffraction region twice, with the same value of an- 
gular momentum. The probability to find two points of 
trajectory of length n in the certain part of the phase 
space growth like n 2 and this is expected behavior of 
all terms in the right hand side of Eqs. (^). Because of 
Eq. (|ll|) these terms have different signs and the overall 
result behaves like n. There is a similar effect in theory 
of disordered metals. The interference correction to the 
density evolution operator consists from three evolution 
operators connected by the kernel. This kernel is called 
the Hikami boxO; it is proportional to the small factor ujt 
because of the well known cancelation!!! 
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II. LEVEL STATISTICS 



T > T* , 



(17) 



The form-factor of the two-point correlation function 
of enepgy levels is expressed in terms of scattering matrix 
tracesEJ 



K(r) = ^\tvS 

T = Tl/Nh 



n|2 



n 

— tr PD" 
N H 1 



fl2) 



In the first order in r one obtains K(t) = 2r, because 
of Eq. (Q)L_.this universal result is a consequence of the 
ergodicityEl 

We will go beyond this approximation by taking into 
account the interference correction to the both diffusion 
and cooperon evolution operators. We should introduce 
the third kind of evolution operators, which would take 
into account exactly the diffraction corrections Eq. (||). 
Let us, instead, include these corrections into the defini- 
tion of the operators D and C . Then the next order in r 
correction to the form-factor becomes 

8K(r)=rn J d3 dA £ C^o (3; 4) (13) 
x (^'(4;3)2ReF / (3;4) + ££=o' (3; 3) F D {3; 4) 



The first term in the braces is the sum over periodic or- 
bits of length n visiting twice the diffraction region with 
the same value of angular momentum. The probability 
to find such a trajectory is 



n(n — 1) a 1 

2 fiNH ' 



(14) 



where a = J dld2F D (l,2) w (kR) 1/3 is the effective 
phase volume of the diffraction region, and the f2 is the 
volume of the phase space; it is just twice the perimeter 
of the billiard il — 8ka. 

The second term in the braces in Eq. ( |i"3| ) is the sum 
over periodic orbits which has a self-tracing piece. The 
probability to find such a trajectory is 



n(n — 1) a 



rj 1 



-ja/Q 



(15) 



where j is the length of the self-tracing part of the tra- 
jectory. The first term in the braces in Eq. (|l3|) partially 



cancels the second term, because of Eq. (|11|). Averaging 
over j and assuming n> 1 we obtain 



K(t) = 2t - r z 1 



0(r 3 



(16) 



where the total amount of periodic trajectories of the 
length n was calculated from the sum rule Eq. (^). The 
probabilistic treatment of the right hand side of Eq. ( p"3| ) 
is justified for large enough n, when only the zero mode 
of the classical evolution operators is important. This 
condition can be written as 



where = n^/Nn is the time of mixing being measured 
in the units of the Heisenberg time. At the moment it is 
not clear whether we should take in to account the corre- 
lation between the classical trajectories on the Ehrenfest 
time scalea. 

The linear in r term in the right hand side of 
Eq. (|l6|) implies level repulsion and it cannot be cor- 
rect if R is small and one computes the energy levels 
perturbativelyn. The perturbation theory works when 
kR < (9 • 2 6 tt 2 ) 1 / 8 , and therefore Eq. © is valid un- 
der the condition kR ^ 1. This implies a S> 1. Since 
a oc ft -1 / 3 and r» oc ft, our result Eq. ( |l6|) demonstrates 
two regimes 

K(t) = It - ^-ar 3 + 0((ar) 2 ) r, <C r < a' 1 , (18) 



2il 
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K{t) = 2t - t 2 + 0(t 



a < t < 1 



(19) 



The last expression matches the universal form-factor 
K{r) = 2r - r 2 + r 3 /2 + 0(r 4 ). The universality 
of the forrtufactor of classically chaotic systems was 
conjectured! 3 ] on basis of numerical data. Our theory 
supports this conjecture under the condition of Eq. jis|). 

The interference correction to the form-factor, 
Eq. (fl8|), can be-.obtained on the language of the ac- 
tion correlationail One should pair off the trajectory 
with two small angle scatterings and the nearby trajec- 
tory without these two scatterings. Each small angle re- 
flection changes the phase of the wave by ir/3, then the 
total phase difference between two periodic trajectories 
is 27r/3. This pair contributes e 27 ™/ 3 to the probability 
to return. The time reversal pair contributes e~ 27 ™/ 3 and 
the sum is just —1, that is the right sign of the interfer- 
ence correction. 

The special attention must be paid to the case of 
the mixed boundary conditions at the opening of the 
billiard. £3 The phase of the exact scattering matrix is 



4>i(k) = 2arg Hi{kR) - nH[(kR) 
= (j>i(0) - 2arctanL^(0)/2 



(20) 



where Hi and H[ are the Hankel function and its deriva- 
tive, k is the parameter, the degree of the mixing. Since 
fa —dcj)/dl one get 



w i+K (O) 



(21) 



Therefore the mixed boundary conditions just shift the 
diffraction edge leaving the phase at the tangency un- 
changed. For the large values of k the exact position of 
the diffraction edge is given by complicated expression, 
but the physical mechanism of the interference remains 
essentially the same. 
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III. GENERAL DISCUSSION 



Present work breaks the common believe that the 
level statistic-jof the ergodic systems is determined by 
symmetryJlS'tj This believe is a result of the analyzea-of 
the effective Lagrangian in theory of disordered solidsll3. 
At the wave number equal Uxxero the interaction of the 
diffusion modes is universalEfrEj and implies the univer- 
sality of the level statistics. 

However, this universality is not occasional, it is the 
consequence of the ergodicity and the particle number 
conservation Lajw-n Therefore the interaction of the Liou- 
villian modestiE3 in chaos should be universal too. For 
example one may introduce the ^-correlated disorder po- 
tential as the source of the interference in a chaotic sys- 
tem and obtain the interaction of the._Liouvillian modes 
independent of the potential strength.tJ The same result 
may be obtained in the model with the smooth disorder 
potentials. 

In this work we consider an example of chaotic system, 
the Sinai billiard, in absence of any disorder potential. In 
such a system, the interference between classical trajec- 
tories takes place because of the hard wall diffraction. In 
the semiclassical limit the phase volume of diffraction is 
relatively small. The interaction of the Frobenius-Perron 
modes (we have introduced the area preserving map in- 
stead of the Hamiltonian flow) is proportional to this 
volume and it is not universal. However, the relevant 
parameter is the time Td, a particle needs to enter the 
region of diffraction. At low frequencies utjj < 1, or 
long times f>T/), the density evolution operator is not 
any more Frobenius-Perron, because of diffraction. Ac- 
cording to our results, the interaction between modes of 
this operator is universal. 

The domain of the universality of the interference ef- 
fect is, therefore, A/h <C u -C t^ 1 and it depends ex- 
plicitly on the phase volume of diffraction. In the case 
of the Sinai billiard this domain is large enough and our 
result supports the-.universality of the level statistics ob- 
served numericallyatJ. Generic chaotic system can have 
so long td at given energy, that its level statistics will 
never manifest the universality. 

In summary we have found that the hard wall diffrac- 
tion contributes the r 2 and possibly the high order terms 
into the form-factor of the energy levels correlation func- 
tion of the Sinai billiard. 
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